3.1.95 \(\int x^3 (A+B x^2) (a+b x^2+c x^4)^3 \, dx\) [95]

3.1.95.1 Optimal result
3.1.95.2 Mathematica [A] (verified)
3.1.95.3 Rubi [A] (verified)
3.1.95.4 Maple [A] (verified)
3.1.95.5 Fricas [A] (verification not implemented)
3.1.95.6 Sympy [A] (verification not implemented)
3.1.95.7 Maxima [A] (verification not implemented)
3.1.95.8 Giac [A] (verification not implemented)
3.1.95.9 Mupad [B] (verification not implemented)

3.1.95.1 Optimal result

Integrand size = 25, antiderivative size = 166 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{4} a^3 A x^4+\frac {1}{6} a^2 (3 A b+a B) x^6+\frac {3}{8} a \left (a b B+A \left (b^2+a c\right )\right ) x^8+\frac {1}{10} \left (3 a B \left (b^2+a c\right )+A \left (b^3+6 a b c\right )\right ) x^{10}+\frac {1}{12} \left (b^3 B+3 A b^2 c+6 a b B c+3 a A c^2\right ) x^{12}+\frac {3}{14} c \left (b^2 B+A b c+a B c\right ) x^{14}+\frac {1}{16} c^2 (3 b B+A c) x^{16}+\frac {1}{18} B c^3 x^{18} \]

output
1/4*a^3*A*x^4+1/6*a^2*(3*A*b+B*a)*x^6+3/8*a*(a*b*B+A*(a*c+b^2))*x^8+1/10*( 
3*a*B*(a*c+b^2)+A*(6*a*b*c+b^3))*x^10+1/12*(3*A*a*c^2+3*A*b^2*c+6*B*a*b*c+ 
B*b^3)*x^12+3/14*c*(A*b*c+B*a*c+B*b^2)*x^14+1/16*c^2*(A*c+3*B*b)*x^16+1/18 
*B*c^3*x^18
 
3.1.95.2 Mathematica [A] (verified)

Time = 0.04 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.00 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{4} a^3 A x^4+\frac {1}{6} a^2 (3 A b+a B) x^6+\frac {3}{8} a \left (a b B+A \left (b^2+a c\right )\right ) x^8+\frac {1}{10} \left (3 a B \left (b^2+a c\right )+A \left (b^3+6 a b c\right )\right ) x^{10}+\frac {1}{12} \left (b^3 B+3 A b^2 c+6 a b B c+3 a A c^2\right ) x^{12}+\frac {3}{14} c \left (b^2 B+A b c+a B c\right ) x^{14}+\frac {1}{16} c^2 (3 b B+A c) x^{16}+\frac {1}{18} B c^3 x^{18} \]

input
Integrate[x^3*(A + B*x^2)*(a + b*x^2 + c*x^4)^3,x]
 
output
(a^3*A*x^4)/4 + (a^2*(3*A*b + a*B)*x^6)/6 + (3*a*(a*b*B + A*(b^2 + a*c))*x 
^8)/8 + ((3*a*B*(b^2 + a*c) + A*(b^3 + 6*a*b*c))*x^10)/10 + ((b^3*B + 3*A* 
b^2*c + 6*a*b*B*c + 3*a*A*c^2)*x^12)/12 + (3*c*(b^2*B + A*b*c + a*B*c)*x^1 
4)/14 + (c^2*(3*b*B + A*c)*x^16)/16 + (B*c^3*x^18)/18
 
3.1.95.3 Rubi [A] (verified)

Time = 0.45 (sec) , antiderivative size = 170, normalized size of antiderivative = 1.02, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.120, Rules used = {1578, 1195, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx\)

\(\Big \downarrow \) 1578

\(\displaystyle \frac {1}{2} \int x^2 \left (B x^2+A\right ) \left (c x^4+b x^2+a\right )^3dx^2\)

\(\Big \downarrow \) 1195

\(\displaystyle \frac {1}{2} \int \left (B c^3 x^{16}+c^2 (3 b B+A c) x^{14}+3 c \left (B b^2+A c b+a B c\right ) x^{12}+\left (B b^3+3 A c b^2+6 a B c b+3 a A c^2\right ) x^{10}+\left (3 a B \left (b^2+a c\right )+A \left (b^3+6 a c b\right )\right ) x^8+3 a \left (a b B+A \left (b^2+a c\right )\right ) x^6+a^2 (3 A b+a B) x^4+a^3 A x^2\right )dx^2\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{2} \left (\frac {1}{2} a^3 A x^4+\frac {1}{3} a^2 x^6 (a B+3 A b)+\frac {3}{7} c x^{14} \left (a B c+A b c+b^2 B\right )+\frac {3}{4} a x^8 \left (A \left (a c+b^2\right )+a b B\right )+\frac {1}{6} x^{12} \left (3 a A c^2+6 a b B c+3 A b^2 c+b^3 B\right )+\frac {1}{5} x^{10} \left (A \left (6 a b c+b^3\right )+3 a B \left (a c+b^2\right )\right )+\frac {1}{8} c^2 x^{16} (A c+3 b B)+\frac {1}{9} B c^3 x^{18}\right )\)

input
Int[x^3*(A + B*x^2)*(a + b*x^2 + c*x^4)^3,x]
 
output
((a^3*A*x^4)/2 + (a^2*(3*A*b + a*B)*x^6)/3 + (3*a*(a*b*B + A*(b^2 + a*c))* 
x^8)/4 + ((3*a*B*(b^2 + a*c) + A*(b^3 + 6*a*b*c))*x^10)/5 + ((b^3*B + 3*A* 
b^2*c + 6*a*b*B*c + 3*a*A*c^2)*x^12)/6 + (3*c*(b^2*B + A*b*c + a*B*c)*x^14 
)/7 + (c^2*(3*b*B + A*c)*x^16)/8 + (B*c^3*x^18)/9)/2
 

3.1.95.3.1 Defintions of rubi rules used

rule 1195
Int[((d_.) + (e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))^(n_.)*((a_.) + (b_.)*(x 
_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(d + e*x)^m*(f + 
 g*x)^n*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, n}, x 
] && IGtQ[p, 0]
 

rule 1578
Int[(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_ 
)^4)^(p_.), x_Symbol] :> Simp[1/2   Subst[Int[x^((m - 1)/2)*(d + e*x)^q*(a 
+ b*x + c*x^2)^p, x], x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x] && Int 
egerQ[(m - 1)/2]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.1.95.4 Maple [A] (verified)

Time = 0.15 (sec) , antiderivative size = 170, normalized size of antiderivative = 1.02

method result size
norman \(\frac {a^{3} A \,x^{4}}{4}+\left (\frac {1}{2} A \,a^{2} b +\frac {1}{6} B \,a^{3}\right ) x^{6}+\left (\frac {3}{8} A c \,a^{2}+\frac {3}{8} A a \,b^{2}+\frac {3}{8} B \,a^{2} b \right ) x^{8}+\left (\frac {3}{5} A a b c +\frac {1}{10} A \,b^{3}+\frac {3}{10} a^{2} B c +\frac {3}{10} B a \,b^{2}\right ) x^{10}+\left (\frac {1}{4} A a \,c^{2}+\frac {1}{4} A \,b^{2} c +\frac {1}{2} B a b c +\frac {1}{12} B \,b^{3}\right ) x^{12}+\left (\frac {3}{14} A b \,c^{2}+\frac {3}{14} B a \,c^{2}+\frac {3}{14} B \,b^{2} c \right ) x^{14}+\left (\frac {1}{16} A \,c^{3}+\frac {3}{16} B b \,c^{2}\right ) x^{16}+\frac {B \,c^{3} x^{18}}{18}\) \(170\)
gosper \(\frac {1}{4} a^{3} A \,x^{4}+\frac {1}{2} x^{6} A \,a^{2} b +\frac {1}{6} x^{6} B \,a^{3}+\frac {3}{8} x^{8} A c \,a^{2}+\frac {3}{8} x^{8} A a \,b^{2}+\frac {3}{8} x^{8} B \,a^{2} b +\frac {3}{5} x^{10} A a b c +\frac {1}{10} x^{10} A \,b^{3}+\frac {3}{10} x^{10} a^{2} B c +\frac {3}{10} x^{10} B a \,b^{2}+\frac {1}{4} x^{12} A a \,c^{2}+\frac {1}{4} x^{12} A \,b^{2} c +\frac {1}{2} x^{12} B a b c +\frac {1}{12} x^{12} B \,b^{3}+\frac {3}{14} x^{14} A b \,c^{2}+\frac {3}{14} x^{14} B a \,c^{2}+\frac {3}{14} x^{14} B \,b^{2} c +\frac {1}{16} x^{16} A \,c^{3}+\frac {3}{16} x^{16} B b \,c^{2}+\frac {1}{18} B \,c^{3} x^{18}\) \(194\)
risch \(\frac {1}{4} a^{3} A \,x^{4}+\frac {1}{2} x^{6} A \,a^{2} b +\frac {1}{6} x^{6} B \,a^{3}+\frac {3}{8} x^{8} A c \,a^{2}+\frac {3}{8} x^{8} A a \,b^{2}+\frac {3}{8} x^{8} B \,a^{2} b +\frac {3}{5} x^{10} A a b c +\frac {1}{10} x^{10} A \,b^{3}+\frac {3}{10} x^{10} a^{2} B c +\frac {3}{10} x^{10} B a \,b^{2}+\frac {1}{4} x^{12} A a \,c^{2}+\frac {1}{4} x^{12} A \,b^{2} c +\frac {1}{2} x^{12} B a b c +\frac {1}{12} x^{12} B \,b^{3}+\frac {3}{14} x^{14} A b \,c^{2}+\frac {3}{14} x^{14} B a \,c^{2}+\frac {3}{14} x^{14} B \,b^{2} c +\frac {1}{16} x^{16} A \,c^{3}+\frac {3}{16} x^{16} B b \,c^{2}+\frac {1}{18} B \,c^{3} x^{18}\) \(194\)
parallelrisch \(\frac {1}{4} a^{3} A \,x^{4}+\frac {1}{2} x^{6} A \,a^{2} b +\frac {1}{6} x^{6} B \,a^{3}+\frac {3}{8} x^{8} A c \,a^{2}+\frac {3}{8} x^{8} A a \,b^{2}+\frac {3}{8} x^{8} B \,a^{2} b +\frac {3}{5} x^{10} A a b c +\frac {1}{10} x^{10} A \,b^{3}+\frac {3}{10} x^{10} a^{2} B c +\frac {3}{10} x^{10} B a \,b^{2}+\frac {1}{4} x^{12} A a \,c^{2}+\frac {1}{4} x^{12} A \,b^{2} c +\frac {1}{2} x^{12} B a b c +\frac {1}{12} x^{12} B \,b^{3}+\frac {3}{14} x^{14} A b \,c^{2}+\frac {3}{14} x^{14} B a \,c^{2}+\frac {3}{14} x^{14} B \,b^{2} c +\frac {1}{16} x^{16} A \,c^{3}+\frac {3}{16} x^{16} B b \,c^{2}+\frac {1}{18} B \,c^{3} x^{18}\) \(194\)
default \(\frac {B \,c^{3} x^{18}}{18}+\frac {\left (A \,c^{3}+3 B b \,c^{2}\right ) x^{16}}{16}+\frac {\left (3 A b \,c^{2}+B \left (a \,c^{2}+2 b^{2} c +c \left (2 a c +b^{2}\right )\right )\right ) x^{14}}{14}+\frac {\left (A \left (a \,c^{2}+2 b^{2} c +c \left (2 a c +b^{2}\right )\right )+B \left (4 a b c +b \left (2 a c +b^{2}\right )\right )\right ) x^{12}}{12}+\frac {\left (A \left (4 a b c +b \left (2 a c +b^{2}\right )\right )+B \left (a \left (2 a c +b^{2}\right )+2 b^{2} a +c \,a^{2}\right )\right ) x^{10}}{10}+\frac {\left (A \left (a \left (2 a c +b^{2}\right )+2 b^{2} a +c \,a^{2}\right )+3 B \,a^{2} b \right ) x^{8}}{8}+\frac {\left (3 A \,a^{2} b +B \,a^{3}\right ) x^{6}}{6}+\frac {a^{3} A \,x^{4}}{4}\) \(226\)

input
int(x^3*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x,method=_RETURNVERBOSE)
 
output
1/4*a^3*A*x^4+(1/2*A*a^2*b+1/6*B*a^3)*x^6+(3/8*A*c*a^2+3/8*A*a*b^2+3/8*B*a 
^2*b)*x^8+(3/5*A*a*b*c+1/10*A*b^3+3/10*a^2*B*c+3/10*B*a*b^2)*x^10+(1/4*A*a 
*c^2+1/4*A*b^2*c+1/2*B*a*b*c+1/12*B*b^3)*x^12+(3/14*A*b*c^2+3/14*B*a*c^2+3 
/14*B*b^2*c)*x^14+(1/16*A*c^3+3/16*B*b*c^2)*x^16+1/18*B*c^3*x^18
 
3.1.95.5 Fricas [A] (verification not implemented)

Time = 0.27 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.00 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{18} \, B c^{3} x^{18} + \frac {1}{16} \, {\left (3 \, B b c^{2} + A c^{3}\right )} x^{16} + \frac {3}{14} \, {\left (B b^{2} c + {\left (B a + A b\right )} c^{2}\right )} x^{14} + \frac {1}{12} \, {\left (B b^{3} + 3 \, A a c^{2} + 3 \, {\left (2 \, B a b + A b^{2}\right )} c\right )} x^{12} + \frac {1}{10} \, {\left (3 \, B a b^{2} + A b^{3} + 3 \, {\left (B a^{2} + 2 \, A a b\right )} c\right )} x^{10} + \frac {3}{8} \, {\left (B a^{2} b + A a b^{2} + A a^{2} c\right )} x^{8} + \frac {1}{4} \, A a^{3} x^{4} + \frac {1}{6} \, {\left (B a^{3} + 3 \, A a^{2} b\right )} x^{6} \]

input
integrate(x^3*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x, algorithm="fricas")
 
output
1/18*B*c^3*x^18 + 1/16*(3*B*b*c^2 + A*c^3)*x^16 + 3/14*(B*b^2*c + (B*a + A 
*b)*c^2)*x^14 + 1/12*(B*b^3 + 3*A*a*c^2 + 3*(2*B*a*b + A*b^2)*c)*x^12 + 1/ 
10*(3*B*a*b^2 + A*b^3 + 3*(B*a^2 + 2*A*a*b)*c)*x^10 + 3/8*(B*a^2*b + A*a*b 
^2 + A*a^2*c)*x^8 + 1/4*A*a^3*x^4 + 1/6*(B*a^3 + 3*A*a^2*b)*x^6
 
3.1.95.6 Sympy [A] (verification not implemented)

Time = 0.04 (sec) , antiderivative size = 202, normalized size of antiderivative = 1.22 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {A a^{3} x^{4}}{4} + \frac {B c^{3} x^{18}}{18} + x^{16} \left (\frac {A c^{3}}{16} + \frac {3 B b c^{2}}{16}\right ) + x^{14} \cdot \left (\frac {3 A b c^{2}}{14} + \frac {3 B a c^{2}}{14} + \frac {3 B b^{2} c}{14}\right ) + x^{12} \left (\frac {A a c^{2}}{4} + \frac {A b^{2} c}{4} + \frac {B a b c}{2} + \frac {B b^{3}}{12}\right ) + x^{10} \cdot \left (\frac {3 A a b c}{5} + \frac {A b^{3}}{10} + \frac {3 B a^{2} c}{10} + \frac {3 B a b^{2}}{10}\right ) + x^{8} \cdot \left (\frac {3 A a^{2} c}{8} + \frac {3 A a b^{2}}{8} + \frac {3 B a^{2} b}{8}\right ) + x^{6} \left (\frac {A a^{2} b}{2} + \frac {B a^{3}}{6}\right ) \]

input
integrate(x**3*(B*x**2+A)*(c*x**4+b*x**2+a)**3,x)
 
output
A*a**3*x**4/4 + B*c**3*x**18/18 + x**16*(A*c**3/16 + 3*B*b*c**2/16) + x**1 
4*(3*A*b*c**2/14 + 3*B*a*c**2/14 + 3*B*b**2*c/14) + x**12*(A*a*c**2/4 + A* 
b**2*c/4 + B*a*b*c/2 + B*b**3/12) + x**10*(3*A*a*b*c/5 + A*b**3/10 + 3*B*a 
**2*c/10 + 3*B*a*b**2/10) + x**8*(3*A*a**2*c/8 + 3*A*a*b**2/8 + 3*B*a**2*b 
/8) + x**6*(A*a**2*b/2 + B*a**3/6)
 
3.1.95.7 Maxima [A] (verification not implemented)

Time = 0.19 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.00 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{18} \, B c^{3} x^{18} + \frac {1}{16} \, {\left (3 \, B b c^{2} + A c^{3}\right )} x^{16} + \frac {3}{14} \, {\left (B b^{2} c + {\left (B a + A b\right )} c^{2}\right )} x^{14} + \frac {1}{12} \, {\left (B b^{3} + 3 \, A a c^{2} + 3 \, {\left (2 \, B a b + A b^{2}\right )} c\right )} x^{12} + \frac {1}{10} \, {\left (3 \, B a b^{2} + A b^{3} + 3 \, {\left (B a^{2} + 2 \, A a b\right )} c\right )} x^{10} + \frac {3}{8} \, {\left (B a^{2} b + A a b^{2} + A a^{2} c\right )} x^{8} + \frac {1}{4} \, A a^{3} x^{4} + \frac {1}{6} \, {\left (B a^{3} + 3 \, A a^{2} b\right )} x^{6} \]

input
integrate(x^3*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x, algorithm="maxima")
 
output
1/18*B*c^3*x^18 + 1/16*(3*B*b*c^2 + A*c^3)*x^16 + 3/14*(B*b^2*c + (B*a + A 
*b)*c^2)*x^14 + 1/12*(B*b^3 + 3*A*a*c^2 + 3*(2*B*a*b + A*b^2)*c)*x^12 + 1/ 
10*(3*B*a*b^2 + A*b^3 + 3*(B*a^2 + 2*A*a*b)*c)*x^10 + 3/8*(B*a^2*b + A*a*b 
^2 + A*a^2*c)*x^8 + 1/4*A*a^3*x^4 + 1/6*(B*a^3 + 3*A*a^2*b)*x^6
 
3.1.95.8 Giac [A] (verification not implemented)

Time = 0.27 (sec) , antiderivative size = 193, normalized size of antiderivative = 1.16 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{18} \, B c^{3} x^{18} + \frac {3}{16} \, B b c^{2} x^{16} + \frac {1}{16} \, A c^{3} x^{16} + \frac {3}{14} \, B b^{2} c x^{14} + \frac {3}{14} \, B a c^{2} x^{14} + \frac {3}{14} \, A b c^{2} x^{14} + \frac {1}{12} \, B b^{3} x^{12} + \frac {1}{2} \, B a b c x^{12} + \frac {1}{4} \, A b^{2} c x^{12} + \frac {1}{4} \, A a c^{2} x^{12} + \frac {3}{10} \, B a b^{2} x^{10} + \frac {1}{10} \, A b^{3} x^{10} + \frac {3}{10} \, B a^{2} c x^{10} + \frac {3}{5} \, A a b c x^{10} + \frac {3}{8} \, B a^{2} b x^{8} + \frac {3}{8} \, A a b^{2} x^{8} + \frac {3}{8} \, A a^{2} c x^{8} + \frac {1}{6} \, B a^{3} x^{6} + \frac {1}{2} \, A a^{2} b x^{6} + \frac {1}{4} \, A a^{3} x^{4} \]

input
integrate(x^3*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x, algorithm="giac")
 
output
1/18*B*c^3*x^18 + 3/16*B*b*c^2*x^16 + 1/16*A*c^3*x^16 + 3/14*B*b^2*c*x^14 
+ 3/14*B*a*c^2*x^14 + 3/14*A*b*c^2*x^14 + 1/12*B*b^3*x^12 + 1/2*B*a*b*c*x^ 
12 + 1/4*A*b^2*c*x^12 + 1/4*A*a*c^2*x^12 + 3/10*B*a*b^2*x^10 + 1/10*A*b^3* 
x^10 + 3/10*B*a^2*c*x^10 + 3/5*A*a*b*c*x^10 + 3/8*B*a^2*b*x^8 + 3/8*A*a*b^ 
2*x^8 + 3/8*A*a^2*c*x^8 + 1/6*B*a^3*x^6 + 1/2*A*a^2*b*x^6 + 1/4*A*a^3*x^4
 
3.1.95.9 Mupad [B] (verification not implemented)

Time = 0.04 (sec) , antiderivative size = 169, normalized size of antiderivative = 1.02 \[ \int x^3 \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=x^{10}\,\left (\frac {3\,B\,c\,a^2}{10}+\frac {3\,B\,a\,b^2}{10}+\frac {3\,A\,c\,a\,b}{5}+\frac {A\,b^3}{10}\right )+x^{12}\,\left (\frac {B\,b^3}{12}+\frac {A\,b^2\,c}{4}+\frac {B\,a\,b\,c}{2}+\frac {A\,a\,c^2}{4}\right )+x^6\,\left (\frac {B\,a^3}{6}+\frac {A\,b\,a^2}{2}\right )+x^{16}\,\left (\frac {A\,c^3}{16}+\frac {3\,B\,b\,c^2}{16}\right )+x^8\,\left (\frac {3\,B\,a^2\,b}{8}+\frac {3\,A\,c\,a^2}{8}+\frac {3\,A\,a\,b^2}{8}\right )+x^{14}\,\left (\frac {3\,B\,b^2\,c}{14}+\frac {3\,A\,b\,c^2}{14}+\frac {3\,B\,a\,c^2}{14}\right )+\frac {A\,a^3\,x^4}{4}+\frac {B\,c^3\,x^{18}}{18} \]

input
int(x^3*(A + B*x^2)*(a + b*x^2 + c*x^4)^3,x)
 
output
x^10*((A*b^3)/10 + (3*B*a*b^2)/10 + (3*B*a^2*c)/10 + (3*A*a*b*c)/5) + x^12 
*((B*b^3)/12 + (A*a*c^2)/4 + (A*b^2*c)/4 + (B*a*b*c)/2) + x^6*((B*a^3)/6 + 
 (A*a^2*b)/2) + x^16*((A*c^3)/16 + (3*B*b*c^2)/16) + x^8*((3*A*a*b^2)/8 + 
(3*A*a^2*c)/8 + (3*B*a^2*b)/8) + x^14*((3*A*b*c^2)/14 + (3*B*a*c^2)/14 + ( 
3*B*b^2*c)/14) + (A*a^3*x^4)/4 + (B*c^3*x^18)/18